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Abstract 

We show, by studying in detail the market prices of options on 
liquid markets, that the market has empirically corrected the simple, 
but inadequate Black-Scholes formula to account for two important 
statistical features of asset fluctuations: 'fat tails' and correlations in 
the scale of fluctuations. These aspects, although not included in the 
pricing models, are very precisely reflected in the price fixed by the 
market as a whole. Financial markets thus behave as rather efficient 
adaptive systems. 

Options markets offer an interesting example of the adaptation of a pop- 
ulation (the traders) to a complex environment, through trial and errors and 
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natural selection (inefficient traders disappear quickly). The problem is the 
following: an 'option' is an insurance contract protecting its owner against 
the rise (or fall) of financial assets, such as stocks, currencies, etc. The 
problem of knowing the value of such contracts has become extremely acute 
ever since organized option markets opened twenty five years ago, allowing 
one to buy or sell options much like stocks. Almost simultaneously, Black 
and Scholes (BS) proposed their famous option pricing theory, based on a 
simplified model for stock fluctuations, namely the (geometrical) continuous 
time Brownian motion model. The most important parameter of the model 
is the 'volatility' a, which is the standard deviation of the market price's 
fluctuations. The Black-Scholes model is known to be based on unrealistic 
assumptions but is nevertheless used as a benchmark by market participants. 
Guided by the Black-Scholes theory, but constrained by the fact that 'bad' 
prices lead to arbitrage opportunities, the option market fixes prices which are 
close, but significantly and systematically different from the BS formula. Sur- 
prisingly, a detailed study of the observed market prices clearly shows that, 
despite the lack of an appropriate model, traders have empirically adapted 
to incorporate some subtle information on the real statistics of price changes. 
Although this ability to price financial assets correctly is often assumed in 
the literature (the 'efficient market' hypothesis), it is in general difficult to 
assess quantitatively, because the 'true' value of a stock, if it exists, is diffi- 
cult to determine. The case of option markets is interesting in that respect, 
because the 'true' value of an option is, in principle, calculable. 

More precisely, a 'call' option is such that if the price x(T) of a given 
asset at time T (the 'maturity') exceeds a certain level x s (the 'strike' price), 
the owner of the option receives the difference x(T) — x s . Conversely, if 
x{T) < x s , the contract is lost. To make a long story short |], @, ||, |J, if T is 
small enough (a few months) so that interest rate effects and average returns 
are negligible compared to fluctuations, the 'fair' price C of the option today 
(T = 0), knowing that the price of the asset now is Xq is simply given by M: 

C(X , X a , T) = / dx' (x' - X s )P{x', T\X , 0) (1) 

JXs 

where P(x', T\x , 0) is the conditional probability density that the stock price 
at time T will be equal to x', knowing its present value is x . Eq. (|TJ) means 
that the option price is such that on average, there is no winning party. 
Pricing correctly an option is thus tantamount to having a good model for 



2 



the probability density P(x', T\xo, 0). 

There is fairly strong evidence that beyond a time scale r of the order 
of ten minutes, the fluctuations of prices in liquid markets are uncorrelated, 
but not independent variables |], [7], || |9], |IIJ]. In particular, it has been 
observed that although the signs of successive price movements seem to be 
independent, their magnitude - as represented by the absolute value or square 



of the price increments- is correlated in time p, 10]: this is related to the 



so-called 'volatility clustering' effect [[11], [7|. More precisely one can represent 
the price x(T) of the asset as 

2-1 

x{T) = x + ]T Sx k (2) 

k=0 

where the increments 5xk are obtained as the product of two random variables: 

5x k = e fc 7 fc (3) 

where (ek)k>o is a sequence of independent, identically distributed random 
variables of mean zero and unit variance, and 7^ is a stochastic scale factor 
independent from the e^s. The sequence (7fc)fc>o is considered to be a sta- 
tionary random process but allowed to exhibit non-trivial correlations (see 
below). Under these hypotheses, the conditional distribution of 5xk, condi- 
tioned on 7^, may be written as: 

P(8x k ) ee lp (4) 

where Pq is independent of k. Models with conditionally Gaussian increments 
- i.e. where Pq is a Gaussian - have been extensively studied [[□]] both 
in discrete time (ARCH models) and continuous time (stochastic volatility 
models) settings. The present model is more general since we do not assume 
that Pq is Gaussian. 

Let us first consider the case where 7^ = 70 is independent of k, which 
corresponds to the classical problem of a sum of independent, identically 
distributed variables. Although P(Sx) is strongly non Gaussian (see, e.g. 



[|l"3|j ), it has a finite variance [[| and the Central Limit Theorem |14| tells us 
that for large N — T/r, P(x', T\xq, 0) will be close to a Gaussian. Using 
then Eq. ([![) essentially leads back to the BS formula [0. For finite N, 
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however, there are corrections to the Gaussian, and thus corrections to the BS 
price. More precisely, the difference between P(x', T|xo, 0) and the limiting 
Gaussian distribution G X0t <j 2 may be calculated using a cumulant expansion 



14j. To a very good approximation, the distribution Pq(Sx) is symmetric |L3|, 
U for time scales less than a month i.e. drift effects are negligible compared 
to fluctuations. This in turn implies that the third cumulant, which measures 
the skewness of the distribution, is small, in which case the leading correction 
in the cumulant expansion mentioned above is, for large N, proportional 



to the kurtosis k, defined as k = (5x 4 ) / (Sx 2 ) 2 — 3 [Q. k vanishes if the 
increments 5x are Gaussian random variables, and measures the 'fatness' of 
the tails of the distribution as compared to a Gaussian. 

Neglecting higher order cumulants, the expansion takes the following 
form: 



{p(x',T\x o ,0)-g xOi(7l (x r )}dx' 



-u 2 /2 



(5) 



where u 



x 



x ) 2 /a 2 



Ti 



a\ and being the variance and kurtosis 
is the gaussian centered at xq of variance 



corresponding to the scale T. Q xo ^ 

It is then easy to show, using Eq. (|IJ), that the leading correction to the 
BS price can be reproduced by using the BS formula, but with a modified 
value for the volatility o = <J (5x 2 ) (which traders call the 'implied volatility' 
X), which depends both on the strike price x s and on the maturity T through: 



T,(x s ,T) = a 




- 1 



(6) 



The fact that implied volatility depends on the strike price x s is known as the 
'smile effect', because the plot of S versus x s , for a given value of T = Nr 
has the shape of a smile (see Fig 1). 

That the volatility had to be smiled up was realized long ago by traders 
- this reflects the well known fact that the elementary increments have fat- 
tailed distributions: large fluctuations occur much more often than for a 
Gaussian random walk. 

As shown in Fig. 2, the smile formula @ reproduces correctly the ob- 
served option prices on the 'Bund' market provided the kurtosis kt in formula 
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Figure 1 
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Figure 1: Example of a smile curve: Implied volatility E(x s ,T) vs distance 
from strike price (x — x s ) for a given T. The data shown correspond to all 
227 transactions of December options on the German Bund future (LIFFE) 
on November 13, 1995. This is a very 'liquid' market, meaning that price 
anomalies are expected to be small, in particular for short maturities T. 
Both call and put options are included with put options transformed into call 
options using the put-call parity 0. Volatilities are expressed as annualized 
standard deviation of price differences. According to Eq. (^) the data should 
fall on a parabola. From a fit of the average curvature of this parabola, we 
extract the 'implied kurtosis' K imp for a given N — — . In this particular case 
we find K imp = 1.92 at N — 144 (9 trading days). 
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Figure 2 
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Figure 2: Plot (in log-log coordinates) of the average implied kurtosis fvi mp 
(determined as in Fig. 1) and of the empirical kurtosis kt (determined 
directly from the historical movements of the Bund contract), as a function 
of the reduced time scale N = T/t, t = 30 minutes. All transactions of 
options on the Bund future from 1993 to 1995 were analyzed along with 5 
minute tick data of the Bund future for the same period. The growth of 
the error bars for the latter quantity comes from the fact that less data is 
available for larger N. Finally, we show for comparison a fit with formula 
(D, with g(£) ~ £~ os , which leads to k t ~ T -0 - 6 (dark line). A fit with an 
exponentially decaying g(£) is however also acceptable (dotted line). 
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(|jD becomes itself T- dependent. The shape of the 'implied' kurtosis «i mp (T) 
as a function of T is given in Fig. 2; K imp (T) is seen to decrease more slowly 
than T _1 . However, if the increments Sx were independent and identically 
distributed (i.e. 7^ = 70), one should observe that k t = k/N. 

Let us then study directly the kurtosis of the distribution of the under- 
lying stock, P(x, T\x , 0), as a function of iV = T/t. In Fig. 2, we have also 
shown kt as a function of N. One can notice that not only Tkt is not con- 
stant (as it should if Sx were identically distributed), but actually kt matches 
quantitatively (at least for N < 200) the evolution of the implied kurtosis 
^imp! (Note that there is no adjustable overall factor.) In other words, the 
price over which traders agree capture rather precisely the anomalous evo- 
lution of kt- A similar agreement has been found on other liquid option 
markets, where bid-ask spreads are sufficiently small to ascertain that the 
quoted prices should indeed be set by a fair game condition. For 'over the 
counter' options, this is likely not to be the case, since a rather high risk 
premium is generally included in the price. 

As we shall show now, the non trivial behaviour of kt is related to the 
fact that the scale of the fluctuations 7^ is itself a time dependent random 



variable [H, EL with rather long range correlations. The random character of 
7^ could come from the fact that 7^ is related to the level of market activity, 
which fluctuates with time. 

We define the correlation function of the scale of fluctuations as: 

m _ (Sx 2 k+e 5x 2 k ) - (Sxp 2 

g(£) is normalized such that g(0) = 1. In this case, one can show that Eq. 
@ holds, with kt given by: 



T 
T 



N p 

k t + 6(k t + 2) ^(1 - -)g(£) 
1=1 iV 



(8) 



where k t is the kurtosis of 5x = x(t + r) — x(t). We have computed from 
historical data the correlation function g(£), which we show in Fig. 3. Inter- 
estingly, g[t) decreases rather slowly, as £~ A , with A ~ 0.6±0.1, from minutes 
to several days. A similar decay of g(£) was observed on other markets as 
well, with rather close values for A, such as the S&P500 (for which A ~ 0.37) 
and the dem/$ market (for which A ~ 0.57). Remarkably, Eq. (Rf) with 
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Fig 3 
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Figure 3: Plot of the daily volatility correlation function g(N) for the Bund 
future market, from 1991 to 1995. A fit by g(N) ~ iV~ A with A = 0.6 is 
shown for comparison. The same behaviour is found to persist for intra-day 
fluctuations (see Inset). 

g(£) oc £~ os leads to ht oc T~ os , in good agreement with both the direct de- 
termination of kt and the one deduced from the volatility smile, n\ mp . Note 
that the effect of a non zero kurtosis on Black-Scholes prices was previously 
investigated in |jl5|, [16|. However, the relation between k t and K imp , and their 
anomalous T dependence, were not, to our knowledge, previously reported. 

In conclusion, we have shown by studying in detail the market prices of 
options that traders have evolved from the simple, but inadequate BS formula 
to an empirical know-how which encodes two important statistical features 
of asset fluctuations: 'fat tails' (i.e. a rather large kurtosis) and the fact 
that the scale of fluctuations exhibits slowly decaying correlations. These 
features, although not explicitly included in the theoretical pricing models 
used by traders, are nevertheless reflected rather precisely in the price fixed 
by the market as a whole. Option markets offer an interesting ground where 
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'theoretical' and 'experimental' prices can be systematically compared, and 
were found to agree rather well [17]. This has enabled us to test quantitatively 
the idea that the trader population behaves as an efficient adaptive system. 
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